Spectra of sparse non-Hermitian random matrices: an analytical solution 
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We present the exact analytical expression for the spectrum of a sparse non-Hermitian random 
matrix ensemble, generalizing two standard results in random-matrix theory: this analytical ex- 
pression constitutes a non-Hermitian version of the Kesten-McKay measure as well as a sparse 
realization of Girko's elliptic law. Our exact result opens new perspectives in the study of several 
physical problems modelled on sparse random graphs which are locally tree like. In this context, we 
show analytically that the convergence rate of a transport process on a very sparse graph depends 
in a nonmonotonic way upon the degree of symmetry of the graph edges. 
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Introduction Random matrices are an indispensable 
tool in the study of various problems throughout physics 
and mathematics. They pervade diverse areas as nuclear 
physics number theory 0, quantum chaos and dis- 
ordered mesoscopic systems [3|, high-dimensional statis- 
tics [3l and information theory (sjl, etc. Owing to the 
ubiquitousness of random-matrix inspired models, exact 
analytical results lead to further insights in several fields. 

The lions' share of the research on random-matrix the- 
ory has been dedicated either to dense ensembles or to 
sparse Hermitian ensembles (the latter in the context of 
spectral graph theory). From the abundance of results 
known for dense ensembles we mention here the 

celebrated Wigner's law Q and Girko's elliptic law [lo| 
for, respectively, Hermitian and non-Hermitian matrices. 
A large number of results for spectra of graphs have also 
been derived Probably, the most simple of these 

results is the Kesten-McKay measure for the spectra of 



undirected regular graphs Exact analytical expres- 



sions for spectra of graphs are important for processes on 
graphs: the computation of effective resistances [l^l, syn- 
chronization in the presence of noise mixing times 
of transport processes [isj . bounds on learnin g in infor- 
mation theory [31, models for quan tum chaos 17 1, com- 
binatorial problems on graphs [l8|, thermodynamics of 
crystalline lattices [l^, etc. 

While many exact analytical results for sparse Hermi- 
tian ensembles are known, exact expressions for the spec- 
tra of their non-Hermitian counterparts have not yet been 
derived. This might be due to the difficulty of applying 
the standard manipulations from Hermitian random ma- 
trix theory to the non-Hermitian case 0, @| • We refer to 
[20l . [21I 1 as one of the few works which have develo ped ex- 
act results on sparse non-Hermitian matrices. In i^, 21 1 
Girko's law has been derived for sparse ensembles at high 
connectivities. Despite these efforts, analytical results 
are sparse and at this moment even an equivalent of the 
paradigmatic Kesten-McKay law is not known for non- 
Hermitian ensembles. A non-Hermitian version of the 



Kesten-McKay measure would be an important tool in 
the development of an exact description of certain physi- 
cal processes taking place on graphs with oriented edges. 
This is especially true if this measure would interpolate 
continuously between a fully undirected and a fully di- 
rected graph. As examples of areas which could benefit 
from such a result, we mention biased diffusi on p rocesses 
[i^ l , synchronization , neural networks [24L [isl ] and 
non-Hermitian quantum systems, such as tight-binding 
models with imaginary vector potentials [26l |27| and 
quantum dissipative systems (28j . 

In this work we present an exact analytical formula for 
the spectrum of a sparse non-Hcrmitian random matrix 
ensemble, opening perspectives in several fields which 
benefit from random matrix theory. Such an exact re- 
sult is possible thanks to recent advances in the theory 
of s par se random matrices with a local-tree like structure 



2l|,[3l|,|33|- Our expression reduces to Girko's elliptic law 
in its highly connected limit and to the Kesten-McKay 
law in its Hermitian limit. Hence, it is a non-Hermitian 
equivalent of the Kesten-McKay law and a sparse realiza- 
tion of Girko's elliptic law. Our result is remarkably sim- 
ple and allows to study the evolution of the spectrum as 
a function of the degree of symmetry in the graph edges. 
We illustrate the interest of this non-Hermitian Kesten- 
McKay law through one physical application: the exact 
calculation of the convergence rate of a stochastic diffu- 
sion on a partially-oriented random graph. This process 
can be seen as a toy model for vehicular traffic as well as 
biological transport [j^ . 

The resolvent equations The complex eigenvalues 
{Ai, A2, • • • , Xn} of a non-Hermitian random matrix 
drawn from an ensemble are defined as the roots of the 
polynomial p{X) ~ dot [A^ — XI] of degree TV in A. We 
define the spectrum of the ensemble of random matrices 
An at a, certain point X = x -\- iy hy 



N 



p(X) = lim iV-i V(5 



i=l 



A,(A«)' 



(1) 



2 



assuming self-averaging of this quantity for N oo. 
The spectrum can be formahy related to the resolvent 
G^(A) of A, defined by G^(A) = (A - through 
the equation p{X) = hmjv^oo(-/V7r)~^9*TrG^(A), where 

d* — a-iid (...)* denotes complex conju- 

gation. The above equation follows from the identity 
i9*A^^ = 7TS{x)S{y), which can be proved by integrat- 
ing both sides over a small square centered at the origin. 
When there is no ambiguity, we leave out the subindex 
N in matrices such as A^r. 

The resolvent Gj^{X) is not properly defined at the 
eigenvalues of A, which are distributed over the complex 
plane. Therefore, one cannot apply the usual resolvent 
manipulations, as done for Hermitian matrices (25j . to 
the non-Hermitian case. We can overcome this problem 
through an Hermitization method [1, 0, [s^l- In this 
method one considers the resolvent (ry) of the 2 x 
2N Hermitian block matrix B2N 

= ( A^^- A* ^ 0^ ^ ) ' 

where O^v is an x matrix filled with zeros and r/ e R 
is a regularizer which keeps G^{'q) properly defined on 
the whole complex plane. The resolvent G^(A) follows 
simply from the NxN lower- left block of lim^_^o (77) , 
where G^{rf) are now computed using standard tech- 
niques for Hermitian matrices (3l| . 

We associate a graph to the matrix A^r by con- 
necting all sites with a nonzero Aij element. If 
this graph has a local tree structure, we have the 
following exact expression for the spectrum p{\) = 

-limAr^oo,,,^o(7riV)"^i9*EjIi[G-j]2i, where the 2x2 
matrices {G^} fulfill the closed set of equations, 

G-i = -A(77) - ^ AuGfAu, (3) 
(Gp"V'-Gri+A,,GfA,,, (4) 

withA(,)^(;^ ^JandA,..(;. ) . From 

the point of view of random graphs, the matrix element 
A^i is the weight of the directed edge from node i io i, 
while di contains the indices of the vertices belonging to 
the neighborhood of node i. We could derive Eqs. ^ and 
(1^ from using recursively the Schur-complement formula 
[sH and we refer the reader to [s^l for technical details. 

Equations ^ and @ have been derived for the first 
time by Rogers and Castillo [U using an approach com- 
ing from the theory of spin glasses. For most random 
matrix ensembles, these equations have a highly intri- 
cate structure. To solve them analytically is unfeasible 
without further simplifications. The authors of [2l| have 
confirmed exactness of Eqs. ([5]|3]) by finding an excellent 
agreement with direct diagonalization methods for sparse 
matrices with a local tree structure. In contrast, here 



we provide an analytical solution to these equations for 
a particular graph structure, which leads to an explicit 
formula for p(A). 

In order to extract an interesting solution from Eqs. ([3]) 
and (|4]), capable to interpolate smoothly between a fully 
undirected and a fully directed graph, we now consider 
an ensemble of random sparse non-Hermitian matrices 
Am of dimension NxN having the following constraints 
on its nondiagonal complex entries Aij 

• {A,,,A,^ G {(^+, (0, 0)}, 

with fc > 1 and the diagonal elements An = 0. This 
random ensemble can be represented by a graph of ver- 
tices [l..iV] and directed segments connecting adjacent 
vertices when Ay ^ 0. In the corresponding graph 
every pair of adjacent vertices is connected by two edges 
in such a way that each vertex contains k incoming (out- 
going) edges of weight A+ and k outgoing (incom- 
ing) edges of weight A- (A+). Hence, the disorder in this 
ensemble is due to the topology of the underlying graph, 
while the elements Aij are fixed to cither of two values 
A± . When the entries are real we have a polarized Bethe 
lattice, illustrated in Fig. [T] For A^ = A^^ we recover 
the unorientcd Bethe lattice [35| of degree 2fc, while for 
A_ = we obtain a fully oriented Bethe lattice of in- 
and out-degree equal to k. 




Bethe lattice Polarized Bethe lattice 

(unoriented regular graph) (partially oriented regular graph) 

FIG. 1. (color online). An example of a Bethe lattice and a 
polarized Bethe lattice, both with degree k = 2. The matrix 
elements corresponding with the red (blue) lines are given by 
A+{A.). 

The resolvent equations (O and ^ for this non- 
Hermitian ensemble become 

G^i = -A(77)-fc(A_G+A+ + A+G_A_), (5) 
G±i = G~i + A±G^A^, (6) 

where A± = ^ ^* ^ , and the spectrum follows 

from p(A) = — TT^^ lim,,_j.o c^* [GJji- The simplified equa- 
tions Eqs. © and ^ follow from Eqs. ^ and (g]) 
after considering the transitive structure of the polarized 
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Bcthc ensembles for N —i' oo. Equations ([5]) and ([6|) de- 
termine, for N — >■ oo, the spectrum of a typical matrix 
drawn from our ensemble, since the resolvent equations 
constitute the typical local neighborhoods of our ensem- 
ble. We now determine the solution to the equations (O 
and ©. 

The analytical expression for the spectrum In the 
following we use polar coordinates and set A± = 
P±exp{i9±). Solving Eqs. ([5]) and ([6]) we find the an- 
alytical expressions for G and G±, see supplemental ma- 
terial The spectrum poW for real matrices Ajv 
(with 9+ — 0^ ~ 0) is given by 



2kHp+p_ 



CW 



it} 



c 



(7) 



Q+Q- 



for 



(8) 



and po(A) = otherwise. In Eqs. ([7]) and ([8]) we have de- 
fined the constants H, C, W, S± and Q±, which depend 
upon k and p± as follows: 



2H = 



HpI +p'^) + ^k^pl-pir + 4{k~i)Hp+p.f, 



{pl+pi)H-2{p+p^y 



C = k^{k-l)-^H-^ 

W = [H^ + {2k - l){p+p^ f] , 
Q±=H±{2k-l)p+p_, 

Eqs. (O and ([8]) are the main result of our work. The 
support ([HI) follows from a stability analysis of the trivial 
solution to Eqs. ^ and ((6)). Indeed, for large values of A 
we find one stable trivial solution with po{X) = 0. This 
solution is unstable at the support of po{^), see |3J]. In 



Figs, m and [3] we compare direct diagonalization results 
with, respectively, the analytical expressions for the sup- 
port dSl) and the spectrum ([7]). We find a very good cor- 
respondence in both cases. Similar to Girko's law, the 
support forms an ellipse, but po{X) is non-uniform. 

Remarkably, for complex entries A± {d± ^ 0) the spec- 
trum p(X) and its support follow simply from Eqs. ([7]) 
and ([8]) through a clockwise rotation by an angle 9 = 
(6*+ + 6^)/2 in the (x, 2/)-coordinate system, such that 
p (A) = po (Ae~*^) . Results ioi 6 are visualized in 
Fig. EI 

We discuss below a couple of interesting limiting cases: 
1. Fully directed Bethe lattice: This limit is obtained 
by setting p_ = in Eqs. dH) and ^ 



Pdb(A) = 



k-l 



kp^ 



^ \\\\^-Bpl 



(9) 








FIG. 2. (color online). Direct diagonalization results of ma- 
trices of size 0(16 -I- 3) (symbols) compared with the support 
^ (lines) for fc = 3, p+ = 1 and different values of p_. 




TpfflffflOnfV^ 

X 

FIG. 3. (color online). Three cuts of po(A), Eq. (O, along 
the real direction are compared with direct diagonalization 
results of matrices of size 0{le + 3), averaged over 0(16 -f 4) 
samples. The parameters are: A: = 3, p+ = 1, p_ = 0.5. 



with the support |AP < kp\. This formula has been con- 
jectured before in [3, [s^, but its rigorous proof remains 
an open problem [8|. 

2. Undirected Bethe lattice: By taking the limit p^ — > 
p_ = jp in Eqs. ([7]) and ([8|), we obtain the Kcstcn-McKay 
law [12| for a graph with degree 2k 



Pkm(A) = 5[y) 



k v/V(2fc- 1) 



(10) 



with support |a::| < 2\p\\^2k — 1. For 6* 7^ and p+ P-, 
the supports of the spectra in Fig. [2] reduce to straight 
lines along the major axes of the ellipses. The projected 
spectrum is then given by the Kesten-McKay measure. 
This concentration of all eigenvalues on the straight line 
y = X tan(0) for a given is a general property of matri- 
ces of the type A + exp(z6')A^ € l^)- Therefore, 
the eigenvalues can be brought to the real axis by a rota- 
tion. Non-Hermitian matrices with a real spectrum have 
attracted considerable attention as alternative theories 
to quantum mechanics (STj . 



FIG. 4. (color online). Direct diagonalization results of ma- 
trices of size 0(le + 3) (symbols) compared with the rotated 
support of (|8]) (solid lines) for k = 3, p+ = 1, 0+ = 0, p- — 0.2 
and different values of The orange circle denotes the 

eigenvalues isolated from the bulk, given by k{A+ + A-). We 
notice the non-trivial rotation of the bulk spectrum of an an- 
gle 9, different from the rotation of the isolated eigenvalue of 
an angle 29 around k. 



3. Dense matrices: By rescaling p± — p^/^/k in 
Eqs. d?]) and ([S]), we obtain the highly connected limit 

when fc — > oo 

P(A)^if^)', (11) 



with support 
12 \ _ 



9 ? 
TT \pi^_ — p_ 



7 — ^ — 72 < i ] i ■ When we 



set E (Afj) = 1 and E [^ijAjj] = r, with E(. . . ) denoting 
the ensemble average and —1 < t < 1, Eq. pT|) reduces 
to p(A) = [7r(l-r2)]"^ with support .x^/ (1 + r)^ + 

2/V(1-t)^ < 2 [25|. Thus, Eqs. ^ and © form a 
non-trivial sparse realization of Girko's elliptic law. In- 
deed, for fc — > oo we have highly connected sparse ensem- 
ble leading to Girko's elliptic law. Sparse realizations of 
Wigner's semicircular law and Girko's circular law has 
recently been proven [H, |3^ . 

Biased diffusion on a regular graph As an applica- 
tion we determine the convergence rate to the stationary 
state of a non-equilibrium transport process on a regular 
graph. Consider a set of random walkers moving along 
the edges of our partially-oriented regular graph model 
with 9± ^ Q and transition rates p± > 0. The relative 
occupancies tt = (7ri,7r2,-- - tT^n), with tt^ the relative 
occupancy of the i-th site, fulfill the linear equation 

^TT = Lattt , Ln = An - k{p+ + p^)1n ■ (12) 
at 

Due to the Perron- Frobenius theorem, the Laplacian ma- 
trix L TV has a unique eigenvector with positive entries and 
eigenvalue Aq = 0, which corresponds to the stationary 
solution of Eq. (|12p . All other eigenvalues have a neg- 
ative real part. If we order them as Ao > 72,e(Ai) > 



FIG. 5. (color online). Spectral gap as a function of the 
degree of symmetry a — p-/p+ for p+ = 1 and k given. 
When k — )> oo, the curves become straight lines. 



••• > TZelXpf^i), then we define the spectral gap as 
g = |72,e(Ai)|. This diffusion process converges expo- 
nentially to the steady state at a rate g. 

We use Eq. ^ for the boundary to obtain the spectral 
gap ioT N oo 



9 [h\a,k) + {2k-l)a\ 

— = 1 + a)fc — 

p+ ri[a, fc) 



(13) 



assuming that finite-size effects at the boundaries of 
/9o(A) are negligible for large N Q. The quantity a = 
p^/p^ is the degree of symmetry of the diffusion process 
and 



2h'^{a, k) = fc(l + a^) + v/A:2(l - ^2)2 + ^^j^ _ i^2^2 

(14) 

We plot g as a function of a in Fig. [5] Remarkably, the 
speed of convergence is non-monotonic as a function of 
a for fc = 2, indicating that for intermediate values of a 
the diffusion process converges more slowly to the steady 
state. For fc > 2 the fully symmetric process has the 
fastest convergence. For dense random matrices such as 
the Ginibre ensemble the spectral gap is zero. Hence, we 
have determined analytically a physical property absent 
in previous studies. 

Conclusions We have presented the analytical expres- 
sion for the spectrum of a sparse non-Hermitian random 
matrix ensemble, which reduces to the Kesten-McKay 
law for Hermitian matrices and to the Girko's elliptic 
law in its highly connected limit. Previous studies for 
sparse non-Hermitian random matrix ensembles relied on 
numerical results. Using a specific distribution for the 
matrix elements we have found to our knowledge a first 
analytical expression for a non-Hermitian random matrix 
ensemble generalizing the Kesten-Mckay measure. Such a 
result can stimulate research in two different directions. 
On the one hand, we have presented a random graph 
model which allows to address analytically how various 
processes on networks depend on the orientation of the 
edges. We have illustrated this through a study of the 
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convergence rate of a biased transport process. Exam- 
ples that deserve further study include network synchro- 
nization (23j . neural networks 24, 25| and non-Hermitian 



quantum mechanics j26l - l28j . On the other hand, since 



our result is derived through a heuristic approach, it 
poses a challenge for the development of new rigorous 
methods which deal with sparse non-Hermitian random 
matrices 0]. Further interesting research directions are: 
the development of exact results for eigenvector localiza- 
tion [i^l and distribution of the largest eigenvalue 41] of 
sparse non-Hcrmitian matrices. 
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Supplemental material 1 

Derivation of the resolvent equations 
I. Neri and F. L. Metz 



INTRODUCTION 



Consider a large matrix An of size with matrix elements Aij = [An]-^ G C. In 
this document we present the essential steps in deriving the following set of resolvent 
equations 

Gr^ = -A(77)-$^A,,Gf (1) 
iG?r' = Gri + A,,G«A,„ (2) 

with A(?7) = (^^\* ~ ^ ^ • '^^^ ^« defined as di = 

{j G [l..A^]|Aji 7^ V Aij 7^ 0}. For more details and the definition of other quantities 
we refer to the main paper. 

The interest in the set of equations ([1112]) is that it determines the spectrum oi An 
through 

1 ^ 

p(A) = - hm ^^d*^[G,h,, (3) 

i=l 

where d* = ^ + ^^j- It is also possible to determine other quantities such as the 

diagonal correlator of eigenvectors but we do not elaborate on this point in this 
paper. 

The equations ([T]) and have been derived for the first time in the recent work 
[2|. However, a rigorous derivation of these equations has not been presented yet. The 
conjecture is that they are exact when N ^ oo and when the underlying graph structure 
is locally tree-like. By the underlying graph we mean the connectivity graph one can 
form when connecting all points which have Aij ^ 0. The exactness of Eqs. ([1]) and 
([2]) is strongly supported by a very good agreement between their numerical solution and 
direct diagonalization methods The concept of a local-tree like structure is properly 
defined in [3| and intuitively it means that most of the vertices have a finite neighborhood, 
where short loops are absent. This is true for randomly drawn graphs when they are 
very large Of course it is true that many real- world graphs are not constructed by 



2 

a random algorithm and contain an intricate loop structure. Nevertheless, it remains 
interesting to consider random systems as a first model for complex systems. 

In this supplemental material we present an alternative derivation of the set of equa- 
tions ([Tll2]). The main difference between our method and the approach presented in 
is that we avoid a mapping on a statistical mechanics problem. Therefore, the assump- 
tions become more intuitive as we avoid the factorization of a complex valued function 
[2I, which plays a role analogous to local marginals in spin models. Another advantage 
of our approach is that it stands much closer to resolvent methods used for rigorous 
spectral calculations on dense matrices Q. 



HERMITIZATION PROCEDURE 



The spectral density of An is obtained from the resolvent G^(A) = (A — A)~^ through 
equation p{X) = \imN^ao{N7i)~^d*TTG j^{X) . The first step of the method consists in 
mapping the calculation of the spectrum of the non-Hermitian matrix on a. resolvent 
calculation of the Hermitian matrix B2N B] , defined as 



B 



2N 







TV 



A 



N 



X* 



A 



N 



(4) 



where Oat is a x matrix filled with zeros and G M is a regulator which keeps all 
quantities properly defined. In this framework, the spectrum follows from 



p(A) = -(Nny^d* limTr 



On On 



(5) 



The Hermitization procedure can also be presented graphically, as can be seen in figure 
1. It corresponds to the translation of an oriented graph to an unoriented bipartite 

W 9 W (children) tf— ■ g— ^ 

Directed graph *■ Bipartite 

Hermitization undirected graph 
procedure 

FIG. 1. (color onhne). Left: A non-Hermitian graph with oriented edges. Right: The undirected bipartite graph 
resulting from the application of the Hermitization procedure to the subgraph on the left. 

graph. Indeed, we can associate to the matrix An a, graph G = {V, E, W), with the set 
of vertices V = [1, N] and the set of edges E = {{i,j)\Aij 7^ 0}. Each edge has a weight 
Wij = Aij, which is denoted by the mapping W : E ^ C. The edge is undirected when 



3 

Wij = w*^, otherwise it is directed. Therefore, the translation B2N corresponds 

to a translation from an oriented graph to a bipartite unoriented graph. 



RECURSIVE APPLICATION OF THE SCHUR-COMPLEMENT 

FORMULA 



From equation (??) we see that to calculate the spectrum, we need to have an expres- 
sion for the diagonal elements of the resolvent G = {ril2N — B2n)~^- Since this involves 
the inverse of a matrix, the following formula is very useful: 



A B 
C D 



f {A-BD-^C)~^ -{A-BD-^C) ^ BD'^ \ 

V -D-^C {A - BD-^Cy^ D-^ + D-^C {A - BD-^C)'^ BD-^ J ' 

( 



(6) 



The above formula is referred to as the Schur- complement formula, and it is a common 
tool in the determination of spectra of dense matrices 0] . Recently the Schur- complement 
formula has also been applied to derive the Kesten-Mckay law [6|. Here we present a 
calculation similar to the one presented in 



Bn 




FIG. 2. A graphical representation of the recursive apphcation of the Schur-complement formula on the bipartite 
matrix i?2jv- We present liere the rooted graph around vertex n, corresponding to the matrix i?„. We can close 



the Schur-complement recursion by setting B 



,{m) 



When relabeling the indices of a matrix we do not change its eigenvalues. We use this 
property and relabel all the vertices such that + N) ^ {2i,2i + l). It becomes useful 
to label each couple (2i, 2i + 1) by i and write B — rjl as a. matrix consisting of 2 x 2 

(IT) A \ 
A ir] j ^^^^^ 



the ij-th block is denoted by Aij 



A 

A* 



We therefore find 
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2N 



-A 
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A12 
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(Af-l)l ^(Ar-1)2 



Nl 
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N2 



Ai 

A2{N-1) A 



-A 

^■Ar(Af-l) 



AiN 

2N 



K7V-1)7V 

-A 



(7) 



where the prime denotes the fact that we have relabeled the indices in B2n- In fact, 
the matrix Sgjy is the relabeled matrix following from the above permutation operation 
on the matrix rjl — B2n- Let us consider now the unoriented graph G' = {V, E', W), 
with V = [l..A^], E' = {{i,j)\Aij 7^ 0} and Wij = 1 for all edges. We perform now a 
first depth search around a certain root vertex n in this unoriented graph and label the 
vertices accordingly. We define the rooted matrix S*-"-* = {B' , n), as the matrix which we 
have created from B' through a permutation according to the first-depth search around 
n. We have therefore for S*-"-*: 
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Ann, O2 ■■■ A, 
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nin 
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^n2n 
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(8) 



We have introduced the degree kn, which denotes the number of matrix elements Ani, 
n) G V, which are different than zero. The matrices B^^^ correspond to a subgraph 
of G' associated to one of the children nodes (or branches) of n. This picture is visualized 
in figure 1. 

We apply now the Schur-complement formula, to find: 
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J 11 



J 21 
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J 12 



J 22 



-A - E -II (^nn,02 ■ ■ ■ O2) (Sg)) {Ann.O 



02)' 

(9) 



We see that to determine the right hand side of equation ([H]), it is necessary to have 
an expression in the upper 2x2 dimensional block of the inverse of the matrix -B^j"-*, 



i. e., we need to know the elements 



Til 



J 11 



Til 



J 12 



Til 



-1 



and 



J 21 



si") 



ni 



. These elements can be determined using again the Schur-complement 



J 22 



formula (hence the recursion). Since -B^"-* is the matrix corresponding to one of the 
subgraphs of G', we have (setting ni = m) 
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(10) 



After applying again the Schur-complement formula, we can close the resultant set of 
equations using B^^ = B^^'°\ for all n ^ do and m G do, with do- Using this approxi- 
mation we find the closed set of equations ([1]I2]). Indeed, setting 
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(11) 



and 
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\ V 
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m 
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(12) 



we recover indeed the equations ([T]|2]). 

We remark that the condition B\^ = B^^'"^ is indeed exact on a tree. Therefore, our 
approximation has some clear intuition. 
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Supplemental material 2 

Solving the resolvent equations 
I. Neri and F. L. Metz 



INTRODUCTION 



In this supplemental material we present how the polarized Bethe lattice leads to a 
non-trivial solution to the resolvent equations. We present the exact analytical solution 
to the resolvent equations in this case. 



RESOLVENT EQUATIONS: GENERAL CASE 

Let us consider a sparse non-Hermitian matrix A, with Aij = [A]-j and i,j = l...N. 
We define the set di = {j\Aij 7^ V Aji 7^ 0}. The resolvent equations in the 2x2 
matrices Gj and G^*'' are given by 

Gri = -A(r7)-5^A,,G«A,„ (1) 

eedi 

iGfY' = Gr' + A.,GfA,, (2) 

with X{ri) = ^ and An = ^ ^ . These equations present an algorithm 

which allows us to determine, among other spectral quantities, the spectrum of the graph 
through 



1 ^ 

P(A) = -^ lim ,i^d*J2[Grhi (3) 

i=l 

where d* = ^ + ^^j- The equations ([1]) and ([2]) are conjectured to be exact when 

— 7- 00 and when the graph is locally tree like. Its use as an efficient and accurate 
algorithm has been demonstrated in [l[ . Here we present an analytical solution to Eqs. (P- 
E]) for a partially-oriented regular graph. Our formula generalizes Girko's elliptic law to 
the sparse case as well as the Kesten-Mckay law to the non-Hermitian case. 
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RESOLVENT EQUATIONS FOR PARTIALLY-ORIENTED REGULAR 

GRAPHS 



The resolvent equations for the partially-oriented regular graph, see figure 2 in the 
main paper, are given by 



G-^ = -A(?7) -fc(A_G+A+ + A+G_A_), (4) 
G±i = G-1 + A±G^A^.. (5) 

These are equations (6-7) in the main paper. The set of equations (j4Jl5]) has only one 
stable solution. This has been proven for the Hermitian case [2|, and we conjecture, 
for now, this to be also true in the non-Hermitian case. Something which is supported 
by numerically solving the resolvent equations. We find for large values of |A| a trivial 
solution (G^,G*) (with p(A) = 0). At small values of |A|, the trivial solution becomes 
unstable in favor of a non-trivial solution (G!]_, G'^) (with p(A) > 0). Below we describe 
the trivial solution (G^, G*^) and the non-trivial solution {G^, G"^). We also determine 
the line in the complex plane where the trivial solution becomes unstable. This line 
corresponds with the boundary of the support of the spectrum. 

Let us first define some constants which we will use throughout: 



u = 


A+A*_ + A*^A_ 


2 


V = 


A*^A+ + A*,A_ 


2 


w = 


A.A*_ - A+Al 


2i 


z = 


A+A^ - A\A*_ 


2i 



We remind that A^ and A_ are the non-zero matrix elements or equivalently {A±,A^) 
are the weights of the edges in the graph. Every vertex is incident to k edges of weight 
{A+,A_) and k edges of weight {A^,A+). It is convenient to use polar coordinates and 
write down A^ = exp(i5+) and A^ = p_ exp(z6'_). 

At last we point out that it will be useful to parametrize the resolvent matrices as 
follow: 



with (a, b, c, d) G C^. 
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THE TRIVIAL SOLUTION: p(A) = 



The trivial solution follows from substitution of equations ([MZ]) in dlHS]) and setting 
a = 6 = 0. The trivial solution is therefore given by 
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Co 

rfo 



d'o 



(8) 
(9) 



with 



Co 
do 



2{2k-l){U -iZ) ' 
-A+ v/A2 -4(2fc- 1) (U + iZ) 
2(2fc- 1) (U + iZ) ' 



(10) 
(11) 



and 



2k -1 



d' 



{k - 1) A* + ^{X*f -A{2k-l){U -iZ) 
2k - 1 

(A; - 1) A + -4{2k-l){U -iZ)' 



The spectrum follows from p{X) = —n ^ ofr'^o ~ 0- 



THE NON-TRIVIAL SOLUTION: p(A) > 



4 



The non-trivial solution follows from substitution of equations ([MZ]) in dUE]) and setting 
a 7^ and b ^ 0: 



cd-- 
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2 / ' 



{2k - 1)U X* + {Z{2k - 1) - H)X 

m - {2k - iy{u^ + z^) ' 

-HX* + {-iZ{2k - 1) + {2k - 1)U) A 

iJ2_(2fc-l)2(f/2 + ^2) ^ 



and for the resolvent: 



-X {-{2k - 1)U + H) + Z{2k - l)yy 



{y {{2k - 1)U + H)- Z{2k - l)x 



2 - {2k - 1)2 (f/2 + Z^) 



H 



h' = 0, 



c' = y {2kZH) + X {-H^ - {2k - 1) {U^ + Z^) + 2kUH) 
-I {y + {2k - 1) (^/^ + Z^) + 2kUH) + X {-2kHZ)) 

d! = M-^ \y {2kZH) + X {-H^ - {2k - 1) [U'^ + Z^) + 2kUH) 
+i {y {H^ + {2k - 1) (^/^ + Z^) + 2A;[///) + a; {-2kHZ)) 



{V{H -V)-W' 

w^ + {H- vy 



2n2 



The denominator is given by 
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i72 _ (2A; - 1)2 + Z2) 



(12) 



with 
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M = 4r 



-X {-{2k - 1)U + H) + Z{2k - l)yy 



{V{H-V) - W 

w'^ + {H- vy 



2^2 



+ [y {2kZH) + X {-H^ - {2k - 1){U^ + Z^) + 2A;f/if)] 
+ [y {H^ + {2k - 1){U^ + Z^) + 2kUH) + x {-2kHZ 



We remind that A = x + i?/. The spectrum of the non-trivial solution follows then from 

, d 



p{X) = -n~'^d' 



p (A) = po {Xe-'') , 



with 9 = 9^ + 6^ and 



2kHp^p^ 
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with constants 



(14) 



(13) 



2H = k{pi + pi) + ^k^{pi - ply + A{k - 1)2 {p+p_)\ 

C = k\k - ir'H-' [{pI + p'_)H - 2 {p^p^f] , 
W= [H' + {2k-l){p^p,y], 
Q± = H±{2k-l)p+p., 
Si = Ql[{HTP-,P-r-HC]- 



I -1 



STABILITY ANALYSIS: DERIVATION OF THE SUPPORT 



The support is found through a stability analysis around the trivial solution (^M)- We 
make therefore an expansion around the trivial solution of the kind 



a = ai, 

C = Co + Ci , 

d = do + di, 



(15) 
(16) 
(17) 
(18) 



with |ai| ^ 1, 1 61 1 <^ 1, |ci| <^ 1 and \di\ <^ 1. Substitution in equations (HHS]), this 
hnear stabihty analysis leads to the following equation: 

\{cl + dl)\ = H-' (19) 

The expression in (??) is in fact an ellipse of the kind: 

+ QZ^ < H~\ (20) 

for ^ = 0, while for ^ > the ellipse rotates around the origin by an angle 6. This 
is the formula presented in our main paper. One can see that Eq. (fT9|) represents an 
ellipse rather quickly by considering that the conformal map C = A± v/F^ maps a 
circle in the (^-space on an ellipse in the A-space. Considering the formula for cq and do, 
Eqs. fllOtllip . we notice indeed that (1191) represents an ellipse. 
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